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ABSTRACT 

We present a semi-analytical model to test the hypothesis that Lya QSO absorption lines 
originate in gaseous halos produced by multiple supernova explosions occurring in Pop III 
objects in a CDM cosmological scenario hoping to assess the validity of CDM models and/or 
constrain their parameters. The preliminary results indicate that the range of Nhi, redshift 
distribution and metallicity of clouds are well reproduced by CDM if they are associated with 
galaxy halos or groups. Firm conclusions on clouds with Nhi < 10 14 cm~ 2 need a more refined 
study. 

Subject headings: Cosmology: theory - dark matter -quasars: absorption lines 



- 3 - 



1. Introduction 

One of the most challenging current attempts in cosmology is to relate the origin and evolution of 
large-scale structure to the observed properties of QSO absorption-line systems. Several groups have now 
performed sophisticated hydro simulations of various hierarchical models (CDM, CDM+A, HCDM) to 
clarify the above issues. In spite of the great advances provided by numerical approaches, their results 
are still not conclusive (and often hard even to compare) since the physics included is necessarily limited 
by computational economy. In particular, the study of effects as feedback and radiative transfer, whose 
importance, particularly for small objects - such as, presumably, the very first generation of galaxies 
(Pop III) - can hardly be overlooked, is still in its youth (Gncdin & Ostriker 1997). These facts motivate 
the effort to develop versatile semi-analytical models, that could co-operatively complement numerical 
studies. Here we present one such whose aim is to test the hypothesis that Lya QSO absorption lines 
originate in gaseous halos produced by multiple supernova explosions occurring in Pop III objects in a 
CDM cosmological scenario. 



2. Dark halos in CDM models 



We assume a CDM power spectrum of fluctuations as in Efstathiou et al. (1992), 
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where ^=1.13; B~6A h~ 2 Mpc; C=3.0 h~ 2 Mpc; D=1.7 hr 2 Mpc; n=l. The previous spectrum describes 
a "standard" CDM model with flo = 1 and h = 0.5, our adopted values; we also take Sit, = 0.19 (from Coma 
cluster data, White & Frenk 1991). The power spectrum normalization has been determined from COBE, 
yielding <r% = 1.07 (Efstathiou et al. 1992). To obtain the number density of dark matter (DM) halos, n, 
as a function of redshift, z, and of their circular velocity, v c = [GM{r) /r] 1 / 2 , we use the Press - Schechter 
formalism (Press & Schechter 1974, PS) following White & Frenk (1991): 
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S c = 1.68 is a critical value of the overdensity and A 2 (r ) is the variance of fluctuations on scale r : 
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-4 - 



We used three different window functions W(kr ). (i) gaussian, (ii) sharp k-space, (iii) top-hat. The results 
are very weakly dependent on the specific choice of W(kr ); we adopt a top-hat window function since it 
gives the best agreement with numerical simulations for 6 C = 1.68 (Lacey & Cole 1994). 



As DM halos arc formed, baryons settle dissipatively into their potential well, being heated by 
virialization shocks. Both baryonic evolution and galaxy formation can be influenced by another important 
heating/ionization source, the ultraviolet background radiation field (UVB); its potential importance has 
been recently discussed by Efstathiou (1992) (see also Thoul & Weinberg 1996), who emphasized its role in 
delaying or even suppressing the formation of very low v c objects. 



In the absence of ionizing radiation, baryons can collapse freely in the gravitational potential of the 
DM, allowing for the formation of the first proto-galaxies. These objects, which we will refer to collectively 
as Pop III, should have masses Mp £ 10 6 M Q in order to be able to cool in a Hubble time (Tcgmark et 
al. 1997). According to PS, the formation peak for these objects is at Zf ~ 30. Assuming an efficiency 
of UV photon production per collapsed proton ~ 190 (Tegmark et al. 1994), and an average life time of 
massive stars ~ 10 7 yr, the ionizing photon rate per galaxy is ~ 7 x 10 50 s _1 . These photons will produce 
expanding HII regions in the surrounding IGM, eventually overlapping. To determine the epoch of complete 
reionization we calculate the IGM porosity parameter, Q, defined by dQ{z) = \dn(Mp, z)/dz\ri(z) 3 dz, with 
Ti ~ (1 + z)~ 3//2 Mpc being the radius of the HII region. Complete overlap (i.e. Q ~ 0.16, Smith 1976) 
occurs at z io „ ~ 7; this value depends on Mp (z ion <~ 4.5[8] for Mp = 5 x 10 6 [5 10 5 ]M Q ). For z < z ion the 



UVB is approximated as a power law spectrum with index a and intensity at the Lyman limit (Haardt & 
Madau 1996, Ferrara & Giallongo 1996 [FG]): 



3. Evolution of the baryonic component 



3.1. Reionization epoch and UVB 







z > 7. 



J_ 21 (z) = 0.53 < 



cxp[-0.69(z-3)] 3 < z < 7, 
(1 + z) 1 / 4 2<z<3, 
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where J_2i = Jo/10 21 ergs s 1 cm 2 Hz 1 sr 
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3.2. Galaxy formation 

As the UVB heats the IGM, it can prevent the gas from collapsing and forming proto-galaxies if the 
IGM temperature, T, is higher than the virial temperature, T c oc^, of the corresponding DM halo. We have 
derived the temperature and ionization evolution of the IGM (assumed of primordial chemical composition) 
from the standard time-dependent energy and ionization equations. The cooling function is described in 
FG; we have studied two cases corresponding to values of a = 1.5(5) mimicking a hard, AGN-dominated 
(soft, galaxy-dominated) spectrum. We have also explored different initial values Tj = T(z = 7) finding 
that this parameter has essentially no influence on the results, as already stressed by FG: we take therefore 
Ti = 10 3 K. The condition T < T c defines a minimum circular velocity, v c ^ m i n {z), (or equivalently a 
minimum mass) for collapse, shown in Fig. 1. Galaxies with total mass below 1 — 4 x 1O 7 M are prevented 
from collapsing after the onset of the UVB; of (M P ) ~ 1.3(M P /10 6 Mo) 1 / 3 km s" 1 

for z > 7. 

3.3. Blow away and secondary halos 

As the gas collapses towards the center of the proto-galaxy, the density enhancement is very likely to 
trigger star formation (SF). Star forming regions will inject energy into the galactic ISM, largely in the 
form of supernova explosions: for a Miller-Scalo IMF, one supernova is expected for each 56 M Q = v~ x 
of stars formed. Type II SNe tend to occur in associations, generating very energetic multiple explosions 
whose observational signature is an expanding supershell. If the galaxy mass is low, its entire gas content 
could be swept by the shell, and driven into the surrounding IGM, in a process called blow-away. The 
expanding shell will be eventually confined by the IGM pressure, provided it is high enough. The blow-away 
condition essentially requires that v s > A(Qo,^b, g) v e, where g is a geometrical factor (details in Ferrara & 
Tolstoy 1997). The previous condition determines the maximum circular velocity for which the blow-away 
can occur, v c>max = 16(1 + z) 1 / 2 km s _1 for our parameters, also shown in Fig. 1; more massive objects will 
escape disruption. Hence, only objects with v c . m i n < v c < v c . max will be able to generate a secondary halo, 
via a shell re-expansion, as explained below. Since the gas is lost, further SF is inhibited, and the remnant 
is a faint agglomerate of low mass stars produced in the first (and last) SF episode. 

Next we derive the evolution of the (spherical) supershell radius, R S} in the thin-shell approximation 
(Kompaneets 1957). The unperturbed external medium has density p e (r) oc r~ 2 for r < r 2 oo, 
corresponding to the galaxy ISM, and p e (r) — p(z) for r > r 2 oo, the IGM density. The shell is 
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driven by the pressure, P b , of the hot gas in the bubble cooling at a rate A(T b ), and subject to the 
gravitational pull of the DM halo (we neglect self-gravity of baryons). The SN mechanical luminosity is 
Lsn oc 7 = vM* ~ {vM gas /Ttf f ) = 1.7 x 10 39 (v c /30 km s" 1 ) 3 ergs s 1 , where 7 is the SN rate, r = 0.6% 
is the SF efficiency and t// is the free-fall time of the gas. The evolution of R s is governed by energy and 
momentum conservation equations (subscripts b and s refer to bubble and shell, respectively): 

^ = L SN - W b P b ^ n 2 b V b A(T b ), (5) 

dt V bpeRs ) = l?7 (n - Pe) - -R^ VbPe > (6) 

where V& = (47r/3)i? 3 . Fig. 2 shows the behavior of R s as a function of z for different values of v c and 
corresponding DM halo formation rcdshift Zf, the evolution is followed up to the stagnation point, i.e. R s ~ 
sound speed in the IGM. Secondary halos associated with DM halos formed at z > 7 are larger, since 
their expansion is unimpeded until the reionization epoch due to the low IGM pressure; the final size 
increases with v c , since L$n oc t> 3 . The final value of R s depends on P e : sizes of the order of tens of kpc 
are obtained setting P e ~ 10 3 x the IGM pressure, or P e /ks ~ 7 K cm~ 3 at z = 0. This value is typical of 
the environments of spiral galaxies and groups: for the Galactic halo at r = 30 kpc, P e /kB ~ 12 K cm~ 3 
(Wolfirc et al. 1995). Lower values of P e result in larger sizes and lower densities. 

Once the stalling point is reached, as the shell pressure is higher than P b , the inner boundary of the 
shell re-expands back into the cavity. The expansion ends when pressure equilibrium across the contact 
discontinuity has been achieved: from this condition we derive the shell thickness, S s . We find that 
the ratio S s /R s > 1/2 in all cases and it often exceeds 0.8. In addition, the re-expansion factor can be 
enhanced by a magnetic field (Slavin & Cox 1992) and by accounting for cooling of hot bubble gas. For 
this reason, we assume that the volume inside R s is uniformly filled with blown-away baryons with density 
n H = m b M/4ir(l - n b )R 3 s . 



3.4. Lya clouds 

We are now able to test the main hypothesis of this work, i.e. that the secondary baryonic halos so 
found can be responsible for Lya absorption lines. As a preliminary test (a complete study is devoted to a 
forthcoming paper) we have derived the neutral hydrogen column density of the clouds and their redshift 
distribution. Assuming a value for the secondary halo gas (or the Lya clouds, in our view) temperature of 
~4x 10 4 K, consistent with a purely thermal Doppler parameter b ~ 25 km s _1 (Lu et al. 1996), we have 
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determined the HI density, nni, from ionization equilibrium, Integrating along a line of sight at impact 
parameter p, we find hydrogen column densities 10 12 < Nhi < 5 x 10 18 cm~ 2 , as shown in Fig. 3, with 
the highest Nhi originating in low Zf halos. Column densities < 10 14 cm -2 are associated with secondary 
halos formed at very high redshift, zj ~ 20, where merging - that we have neglected - should be rather 
important. The redshift distribution of the clouds, 



can be very nicely fit by a power law dN/dz = A (l + z) 1 with 7 = 2.2 and A Q =Q.7 1 in satisfactory 
agreement with observational data (Bechtold 1994, Giallongo et al. 1996) at this stage. 



We have brought arguments suggesting that at least some Lya QSO absorption lines could originate 
in secondary halos produced by multiple supernova explosions occurring in Pop III objects. As a first 
step, we have investigated only a specific CDM model; nevertheless our aim is to extend the analysis to a 
larger class of theoretical models (CDM+A, HCDM) In this paper we have introduced the main features 
of our semi-analytical model, deferring a complete exploration of the results, cosmological models, and 
model details to a forthcoming paper. Yet, two conclusions could be relevant: (i) Lya clouds should be 
preferentially found in relatively high pressure regions, as galaxy halos and groups; (ii) HI column densities 
and redshift distribution are well reproduced by a CDM model if clouds are associated with galaxy halos 
or groups. We reiterate that predictions for Nhi ~ 10 14 cm~ 2 might be affected by the simplifications 
adopted (i.e. no merging). In view of the exploration of additional cosmological models, it is important to 
compare the model predictions on other observables as the dN/dNni and b distributions, and clustering 
properties. The metallicity Z/Z Q ;$ 10~ 2 found in clouds with Nhi > 3 x 10 14 cm~ 2 by Cowie et al. (1995) 
and Tytler et al. (1995) could be explained by primordial enrichment by Pop III objects. A back-on-the 
envelope calculation using the SN rates from our model and assuming a metal yield of ~ 1M Q per SN, gives 
Z/Zq ~ 1.2 x lCT 4 (v c /10 km s" 1 ) 2 in the secondary halo. 





4. 



Conclusions and discussion 
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Fig. 1. — Minimum circular velocity (solid lines) of a DM halo for collapse of the gas as a function of redshift 
for two values of the UVB power-law spectral index a = 1.5,5. Also shown (dashed line) is the maximum 
circular velocity required for blow- away to occur. The initial IGM temperature is Tj = 10 3 K. 



Fig. 2. — Behavior of shell radius, R s , as a function of z for different values of v c (from top to bottom) and 
corresponding DM halo formation redshift zf, the evolution is followed up to the stagnation point. 



Fig. 3. — Nhi of Lya clouds as a function of v c for different values of the DM halo formation redshift Zf. 
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